X-ray tomography of a crumpled plastoelastic thin sheet 
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A three-dimensional X-ray tomography is performed to investigate the internal structure and its 
evolution of a crumpled aluminum foil. The upper and lower bounds of the internal geometric fractal 
dimension are determined, which increase with the compression. Contrary to the simulation results, 
we find that the mass distribution changes from being inhomogeneous to uniform. Corroborated 
with the evidence from previous experiments, these findings support the physical picture that the 
elastic property precedes the plastic one at dominating the deformation and mechanical response 
for all crumpled structures. We show that the interior of a crumpled ball at the plastic regime can 
be mapped to the compact packing of a granular system. 

PACS numbers: 46.32,+x, 62.20.F-, 89.75. Fb, 42.30.Wb 



I. INTRODUCTION 

Crumpling is a common but complicated process that 
occurs in many occasions and at all length scales. It is 
more complex in sheets than wires [l[ because the former 
develops a unique cone structure. Although much has 
been learned from studying the mechanism of a single 
cone[H 01, the collective behavior of vertices is still an 
open field. The main difficulty lies in the treatment of 
topological constraint that they cannot cross each other. 
The first aim of this Letter is to acquire 3D images of 
this exotic self-organized crumpled structure. The con- 
ventional process to achieve this purpose is to unfold and 
profile the topography of its surface, via which the prob- 
ability of different ridge length was found to obey the 
lognormal distribution at large scale and a power law 
at small scale However, so far there is no system- 
atic study on the evolution and statistical properties of 
its labyrinthical interior without destroying the crumpled 
ball. 

Another interesting property is the scaling relation 
proposed by Nelson and KantorQ, who predicted the 
value of dimension D appearing in R oc R where Rq 
and R are the initial size and final radius of the crumpled 
ball. The development of this theory has relied mainly on 
computer simulations over the years. First, Vliegenthart 
and GompperQ studied the elastic sheet and extended 
the scaling relation to include the external force F and 
various modulii: 
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where Kq and k are the two-dimensional Young's mod- 
ulus and bending rigidity, and the exponents a and /3 
are believed to be universal and material-independent. 
More recently, Tallinen et ai.0] included plasticity but 
still found this scaling relation to remain intact. Fur- 
thermore, Astrom et al.&l reported a divergence in F 
when the volume fraction of the sample is around 0.75, 
at which point the elastic behavior turns out to be sim- 



ilar to that of a granular media. In the last decade, sci- 
entists have performed many experiments to determine 
the dimension D of elastic and plastoelastic crumpled 
structures , [ill fl2| . Most were based on the scaling 
relation in Eq.([T]) to give this dimension as a function of 
the two universal exponents, D = 2/(1 + 2/3 — q)(9l. [Tcj| . 
The fractal dimension is believed to be material depen- 
dent but, otherwise, is universal for the external force 
and different initial size, thickness, and final radius of 
the sample. 

It is obvious that a direct observation of the internal 
structure requires a nondestructive method that can see 
through the crumpled ball with sufficient resolution to 
describe the exact location of thin sheets. Normal de- 
structive approach, such as cutting the ball in half to 
expose the cross section is not only affected by possi- 
ble alteration of the detailed structure, but also limited 
severely by the number of planes one can expose and 
therefore reduces the possibility to reconstruct 3D struc- 
tures. X-ray microtomography is certainly the experi- 
mental tool of choice to meet these two demanding re- 
quirements simultaneously. 



II. EXPERIMENTS 

We use the aluminum foils as our sample, which are 
of the same thickness (16/im) but different radius Rq 
[mm] =3, 4, 4.5, 5, 6, 6.5, 7, 8, 9, and 10. They are ran- 
domly folded into a ball of same final radius R — 1.5mm. 
Such a dimension is already difficult to be examined 
by other methods, and hard- X-ray with photon energy 
higher than lOkeV is required. To precisely describe the 
crumpled geometry of the foil with a spacing between ad- 
jacent sheets in the order of /im, high image resolution 
and precision of specimen placement are essential. 

In order not to sacrifice too much resolution, the sam- 
ple is limited to 5 x 5 x 5 mm 3 for the x-ray tomog- 
raphy. This makes the preparation of sample by pres- 
sure chamber fl3l| unrealistic. After trying several tries, 
we eventually settled by mimicking the method first em- 
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ployed by Balankin et al. ; namely, use a flat tip conven- 
tional tweezer to sqeeze while rotating the sample. The 
method is believed to produce reliable and reproducible 
data. 

Under the limitation of beam time, we judge that it 
is more preferable to invest it on samples with diffcr- 
erent parameters. Many anomalous properties have been 
revealed in the literature, which were sensitive to the 
compaction R/Rq of the sample. For instance, there was 
a two-stage transition (folding-crumpling transition) in 
the intermediate range of compaction around 0.23[14j|. 
Furthermore, a jamming transition was found when the 
sample reached the highly compact state. Therefore, a 
detailed scan of different compactions can be informative. 

We employ a special version of X-ray microtomog- 
raphy system based on the high intensity X-ray from 
synchrotron [l5j]. Such systems provide a standard res- 
olution between 1/xm to 2/im which shows clear recon- 
structed images for our analysis. The experiment is per- 
formed at the 01A beamline of National Synchrotron Ra- 
diation Research Center (Hsinchu, Taiwan). The beam- 
line provides unmonochromatic X-rays whose energy dis- 
tribution is 8keV to 15keV. Image acquisition time per 
projection is about 10ms, which is captured by a CCD 
with 2X optical lens focused on a CdWC>4 single crys- 
tal scintillator. The resulted reconstruction consists of a 
data matrix of 1200x1200x1200 pixels with a pixel size 
3/im. Higher resolution image can be obtained by trans- 
mission X-ray micro scop y which has recently achieved 
a resolution of 30nm[16|. However, the high resolution 
comes with a trade-off with reducing specimen size. 




FIG. 1: A slice of raw images reconstructed from 1000 pro- 
jections for R/Ro — 0.167. Before segmentation, two kinds of 
reconstructed artifacts are present. The glisterns are caused 
by the strong refraction and diffraction of incident X-ray from 
aluminum facets when they are aligned to the incident beam 
within a certain value. The ring artifacts are due to the de- 
fects on the detector screen and insufficient background nor- 
malization. Both kinds of noise can be reduced significantly 
by the fill tracing method. 



III. DATA RECONSTRUCTION 

The tomography reconstruction is done with interpo- 
lation on pixels which is not covered during the rota- 
tion. It is therefore important to acquire enough number 
of projections experimentally to perform the reconstruc- 
tion. However, the experimental system, the aluminum 
foil here, is quite simple and its thickness can be assumed 
not to change during the compaction, then the interpo- 
lation is not likely to create problem during the filtered- 
back projection methods for reconstruction. There are 
quite a few more algorithms scientists use to get the to- 
mography reconstruction, but we believe that they do not 
affect the result in this case. The only possible problem 
is that, due to the resolution limit, if the foil wrinkled 
with scale smaller than the resolution then the recon- 
struction will not be able to reproduce it and that could 
affect the numerical analysis. This has been checked by 
opening the foil and confirming that there was no small 
scale distortion. 

In image segmentation, the conventional thresholding 
is an essential and widely-used method to spot the ob- 
ject pixel and convert gray-scale images into binary ones. 
Despite its simplicity, the only tunable parameter, the 
threshold value, constrains the precision of the method. 



Inevitably, some useful information is erased during this 
thresholding of the raw data. To minimize the loss, the 
original images are converted into 1-bit units through the 
fill tracing process (l7j. which makes use of vectorizing to 
trace out the interface between the object and the back- 
ground. To be rid of the tiny areas due to noise and 
dusts, see Fig|TJ minimum area is set at 20 pixels to fil- 
ter out the invalid tracing. Although this method can 
not distinguish two paths when they are in contact, it 
provides a sufficient quality to the information for our 
analysis. All data matrices are reconstructed again af- 
ter the preliminary analysis has been processed. Sample 
pictures are shown in FigO 

The reason why the thresholding is required is that the 
imaging process in the experiment does not give the ab- 
solute gray scale. Every projection is composed of dark 
noises (from CCD, electronic noise, etc.), background 
noises (from the optical elements such as the X-ray win- 
dows, scintillators, lens inferfections, etc.) and other ef- 
fects. Therefore, it is hard to get the absolute density. 
For example, in principle we should be able to extract the 
absorption coefficient of aluminum from the image gray 
scale comparing to air, but that is normally not reliable 
unless a large calibration is done. In this case, we do 
not care about the absolute value and our case becomes 



FIG. 2: (color online) Sequence photos of the reconstructed cross section of crumpled samples. The compactions R/Ro of (a), 
(b), (c) and (d) are 0.38, 0.30, 0.19 and 0.15, respectively. As the compaction decreases, more vertices are induced and more 
locally-aligned structure appears. 



practically a 1-bit system. The only reason to perform 
thresholding is to eliminate noise and artifacts from the 
reconstruction and more subtly the phase contrast effect, 
such as those lines tangential to the edge of the object. 
The danger of thresholding is the elimination of small de- 
tails like those wrinkles. One can answer this concern by 
blowing up the reconstructed image and show the effect 
of different thresholding and its effect on image. 



IV. INTERNAL GEOMETRIC FRACTAL 
DIMENSION 

Fractal dimension is a useful quantity to characterize 
many cluster-assembled structures and the surface of cor- 
rugated thin films with self-affine morphology [HJ]. How- 
ever, cares need to be taken when comparing the inter- 
nal geometric fractal dimension D m with the dimension 
D determined from its mass-size or scaling relation in 
Eq.([T]). For all classic fractals, such that Koch curves, 
Sierpinski gasket and sponge, etc., the equality D m = D 
is obvious from construction. In the case of forced fold- 
ing, the fractal dimension of the set of balls folded from 
sheets of different sizes is determined by the folding con- 
ditions. For example, the set of balls with the same con- 
traction ratio Ro/R has the fractal dimension D = 2. 
At the same time, the sets of balls folded by the same 
force (F= constant) and under the same stress pressure 
(P oc F/R 2 ^constant) are characterized by different 
fractal dimensions Dp > Dp > 2, whereas the inter- 
nal structure of each folded ball does not know how the 
other balls were folded. Furthermore, in the case of (hy- 
per) elastic sheets as, for instance, rubber sheets, the 
balls are completely unfolded after the confinement force 
is withdrawn. In contrast to this, the diameter of folded 
paper sheet increases only slightly increases after with- 
drawing. So, the fractal dimension of the set of folded 
paper balls is dependent on the strain relaxation rate. 

Balankin et qL[Io| have emp loyed the mthod proposed 
by Miyazima and Stanley [lj| to measure D m by pierc- 
ing a needle connected to a string along the diameter of 



the crumpled ball and studying the number of intersec- 
tions as a function of the sheet size. They found that 
the external structures of balls folded from different pa- 
pers are characterized by the same fractal dimension D, 
i.e., D is independent on the mechanical and geometri- 
cal properties of paper sheet. Furthermore, they noted 
that this dimension was numerically close to the value 
of D expected from numerical simulations for the set of 
phantom sheets folded under a fixed force instead of a 
fixed pressure. From this observation they speculated 
that the empirical relation between force and ball diame- 
ter in Eq. ([!} was valid not only for sheets of different size 
Ro, but also for force distribution within the folded ball. 
In their experiments, the ratio Rq/R is relative small in 
the sense that the self-avoiding effects do not play a sig- 
nificant role. When the ball is further compressed, the 
difference between the self-avoiding paper and the phan- 
tom sheets may become more pronouccd. 

Here, we present a more direct measurement of 
D m by adopting the box-counting method to calcu- 
late the Minkowski-Bouligand (MB) dimension, Db = 
-lim £b ^o(logA r fe(eb))/(loge fc ). In FigGH the number of 
occupied boxes, Nb, decays with the box size, e&, in a 
power-law fashion, and the slopes of these solid lines 
are the MB dimension. Since the value of Df, is of- 
ten considered to be the upper bound for the fractal 
dimension, we also calculate the density-density corre- 
lation to determine the correlation dimension to set the 
lower bound. The correlation dimension is defined as 
D c = — lim Cc _ > o(logC(e c )/(loge c ) where the correlation 
function C(e c ) counts the number of pairs (i,j) with 
length s(i,j) < e c . From the right inset in Fig[31 we 
find that both MB and correlation dimensions grow as 
the compaction R/Rq decreases. 



V. MASS DISTRIBUTION 

In the crumpling of an elastic wire, the mass-size rela- 
tion has been studied and compared to the Flory's theory 
for polymers (20l. l2ll. |22| . Adopting Flory's free energy ex- 
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FIG. 3: (color online) The occupied number of boxes is plot- 
ted as a function of box size in a log-log scale. Its slope can 
be readily read off as the MB dimension. Together with the 
correlation dimension, they are shown in the right inset for 
different compaction R/Ro. These two dimensions set the up- 
per and lower bounds for the fractal dimension. We estimate 
it to increase roughly from 2.2 to 2.8 as larger sheets are used 
to achieve lower compactions. The left inset is a log-log plot 
of R versus Ro according to the scaling relation of Eq. {1} . The 
slope of different lines is adopted from the experimental val- 
ues in the right inset. To examine the reproducibility of the 
measurement, a representative error bar in the inset is drawn 
for R/ Rq=Q.Y1 from three more samples. The crumpled ball 
at this ratio just crossed the loose packing regime into the 
compact one. The smallness of error justifies the reliability of 
our data. 



FIG. 4: (color online) The mass density distribution for four 
different compactions R/Ro is plotted as a function of the 
radial length. Inset shows the slope of the main plot. The 
density is divided by Ro rather than R^ for clarity. As the 
compaction decreases, the influence gradually propagates to 
the interior and causes the mass density distribution to ap- 
proach homogeneity. The slope with compaction above 0.43 
is small because the sample has not entered the crumpling 
regime and remains as a sheet. The error bar in the inset is 
drawn for R/ Ro—0.17 from three more samples. 



pression, it was found that the bending and exclusion en- 
ergies dominate the system rather than the entropy (20j. 
This is equivalent to saying that the crumpling is a de- 
terministic process and all the uncertainties of the final 
configuration are attributed to the fluctuation and im- 
perfection during the compression. 

An extensive investigation of the mass distribution 
is done to reveal how the individual mass relaxes lo- 
cally and globally throughout the assembling. A simple 
coarse-grained method is first adopted to calculate the 
mass density distribution. After the radius of gyration, 
^global _ Tmax ^ i s determined, we measure the density 
along the radial direction and average over azimuthal an- 
gles. This simple method, however, suffers from possible 
interference from the deviation of the sample from a per- 
fect sphere. To remedy this defect, we employ another 
more detailed investigation in which the radius of gy- 
ration is defined locally as R l ° c = r l £f ax where r l £f ax is 
determined by the final vanishing of particles within an 
interval of solid angle. We divide the whole sample into 
20 x 36 x 10 pieces along the radial, azimuthal, and tilting 
angles, respectively. After all the masses have been accu- 
mulated in each piece, it is divided by the volume of the 
piece to give the corresponding density, p. In contrast to 
a homogenous distribution, this second method detects a 



linear proportionality for low compaction R/Ro sample, 
see FiglU 

In the inset of FiglU the slope m = dp/dR l ° c decreases 
as R/Rq shrinks. In other words, the mass tends to accu- 
mulate near the crust initially, but eventually shifts ho- 
mogeneously to the core. This trend as well as that of the 
fractal dimension points to the physical picture that elas- 
tic and plastic deformations dominate at separate stages 
during the crumpling. This is contrast to the crumpled 
elastic wires [2(| for which there is no second homoge- 
neous stage. Lin et aL[l3[ have performed 3D crumpling 
experiment and measured the number of layers and ver- 
tices. They found that these two formations dominate at 
different stages. In the beginning, layers are created in 
abundance by the folding. The distribution of facet size 
is wide, which results in a nonuniform mass distribution. 
As the ball size decreases, it becomes difficult to generate 
new layers. Further compression now only serves to bulk 
the existing layers and induce the permanent or plastic 
deformations, such as the ridges and vertices. The small 
facet size that characterizes this stage contributes to a 
more uniform mass distribution. This porous yet com- 
pact packing resembles that of a granular media. This 
analogy has been put forward by Astrom et oi.@| except 
that our transition is gradual and happens much earlier. 
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VI. DISCUSSIONS 

The calculations for both the fractal dimension and 
density exponent involve statistical averaging over the 
bulk of samples. The fact that each crumpled sheet ex- 
hibits thousands upon thousands vertices convinces us 
that it has gone through numerous schochastic processes. 
According to the central limit theorem, the random er- 
rors they create will have cancelled each other consid- 
erably and reached rather accurate properties. Similar 
practice can be found in Ref . [23| . 

Microscopically, we can imagine the thin sheet as a 
network of triangular meshes consisting of springs, such 
as the tethered membrane model by Kantor, Kardar, 
and NelsonQ. In the early stage, the strain is below 
the threshold of plastic deformation and so most springs 
still obey the Hook's stress-strain relation. The correla- 
tion length of the material fully extends and exhibits the 
same property as an elastic sheet. Only when the volume 
fraction is high will the springs enter the plastic regime 
and render the correlation length extremely short. This 
causes the gradual transition into a structure which is 
porous and yet homogenous in mass distribution, similar 
to the granular packing^. 

The volume fraction is estimated to be about 0.45 when 
the compaction R/Rq equals 0.15 which means that more 
than half of the interior is still filled with air. This volume 
fraction is much smaller than the fee packing, 0.74, and 
the resistance divergence point, 0.75, in elastic sheetsQ. 
One question arises: How does this structure derive its in- 
credible resistance from such an inefficient packing com- 
pared to the granular system? One reason is that the 
system is constantly trapped in metastable states 0, [24| 
due to the non-Markovian nature of the process. An- 
other conceivable answer is from the material science. It 
is known in mechanical engineering that a beam buck- 
les under an axial load when its length exceeds 50 to 
200 times of its thickness [25|. For a shorter beam, the 
deformation switches from the bending modulus to the 
normally much-larger Young's modulus. In other words, 
each material has an intrinsic minimum size of facet, be- 
low which it takes an enormously large stress. We believe 
this causes the poor packing of the hollow polyhedrons 
surrounded by these facets inside the crumpled struc- 
ture. The surprising thing is that, rather than compli- 
cating the structure, the plasticity enables the interior to 



reach a higher fractal dimension and more homogeneous 
packing, both of which contribute to its high mechanical 
resistance. 



VII. CONCLUSIONS AND OUTLOOK 

By employing the X-ray microtomography, we find 
that the mass distribution inside the aluminum ball is 
inhomogeneous at low volume fraction and with a fractal 
dimension slightly larger than 2. We call this the loose- 
packing regime. This is reminiscent of the same inhomo- 
geneity observed in the wire crumpling system[l|, 0, [H[ 
which is without the plasticity. As we compress the ball 
further, it enters the compact-packing regime where the 
fractal dimension increases to 2.8 and the mass distribu- 
tion becomes homogeneous. This finding implies that, 
with the shrinkage of ball size, the sheet runs out of 
rooms to generate new layers and thus can only start to 
buckle. When the strain exceeds the yield point, this cre- 
ates a sudden surge of ridges and vertices [13} which mark 
out many separate facets whose size gets smaller as their 
number increases under further compression. The homo- 
geneous mass distribution and strong resistance of this 
configuration are analogous to those of granular packing 
as suggested by previous simulation Q. 

There was a serious debate on whether the scaling law 
holds for real material. Imaging different materials would 
be an exciting outlook. The reason we focused on the 
aluminum foil was due to its stability. Other material, 
such like HDPE, Mylar, and paper, all exhibit a slow 
swelling relaxation [26j. This phenomenon ruins the 3-D 
image reconstruction during the long period of imaging, 
which takes up more than three hours. Adhesive or some 
smart way to confine the sample without introducing fur- 
ther noises may be the future direction to overcome this 
obstacle. 
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